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ABSTRACT 

We investigate the light-element behavior of red giant stars in Northern globular clusters (GCs) 
observed by the SDSS-III Apache Point Observatory Galactic Evolution Experiment (APOGEE). We 
derive abundances of nine elements (Fe, C, N, O, Mg, Al, Si, Ca, and Ti) for 428 red giant stars in 
10 globular clusters. The intrinsic abundance range relative to measurement errors is examined, and 
the well-known C-N and Mg-Al anticorrelations are explored using an extreme-deconvolution code for 
the first time in a consistent way. We find that Mg and Al drive the population membership in most 
clusters, except in M107 and M71, the two most metal-rich clusters in our study, where the grouping 
is most sensitive to N. We also find a diversity in the abundance distributions, with some clusters 
exhibiting clear abundance bimodalities (for example M3 and M53) while others show extended distri¬ 
butions. The spread of Al abundances increases significantly as cluster average metallicity decreases 
as previously found by other works, which we take as evidence that low metallicity, intermediate mass 
AGB polluters were more common in the more metal poor clusters. The statistically significant cor¬ 
relation of [Al/Fe] with [Si/Fe] in M15 suggests that 28 Si leakage has occurred in this cluster. We also 
present C, N and O abundances for stars cooler than 4500 K and examine the behavior of A(C+N+0) 
in each cluster as a function of temperature and [Al/Fe]. The scatter of A(C+N+0) is close to its 
estimated uncertainty in all clusters and independent on stellar temperature. A(C+N+0) exhibits 
small correlations and anticorrelations with [Al/Fe] in M3 and M13, but we cannot be certain about 
these relations given the size of our abundance uncertainties. Star-to-star variations of a—elements 
(Si, Ca, Ti) abundances are comparable to our estimated errors in all clusters. 
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Over the last two decades, the long lasting idea of 
globular clusters hosting single, simple stellar popula¬ 
tions has changed dramatically. The classical paradigm 
of GCs being an excellent example of a simple stellar 
population, defined as a coeval and initially chemically 
homogeneous assembly of stars, has been challenged by 
observational evidence. The presence of chemical in¬ 
homogeneities, in most cases limited to the light ele¬ 
ments (the chemical pairs C-N, O-Na, and Mg-Al anti¬ 
correlated with each other), have been known for decades 
and recognized to be the signature of high-temperature 
H-burning. This was initially fr amed w ithin a stellar evo¬ 
lutionary scenario (see, e.g.. iKraftl I1994L and references 
therein) given that GC abundance work based on high- 
quality data was limited to bright, evolved giants. It was 
only at the turn of the century that the availability of 
high-resolution spectrographs mounted on 8 m class tele¬ 
scopes made it possible to carry out studies on the com¬ 
positions of stars down to the main-sequence, which re¬ 
vealed light ele ments variations an a logous to those found 
among giants (IBrilev et all 119961 IGratton et al.l 120011 : 
IRamirez &; Cohenll2002f) . Given that the atmospheres of 
warm main-sequence stars in large part retain the com¬ 
position of the gas from which they were formed, the 
unavoidable conclusion was that the abundance inhomo¬ 
geneities are of primordial origin. 

The most ext ensive spectroscopic s urvey of GCs un¬ 
dertaken so far (|Carretta et al] l2009a.b.c) revealed that 
these inhomogeneities are ubiquitous in Galactic GCs, 
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though they do not appear to occur in other star for¬ 
mation environments. However, the extent of the in¬ 
homogeneity varies from cluster to cluster, and appears 
to correlate strongly with the present-day total mass of 
the GCs, and also with metallicity. The improvement 
in available instrumentation and techniques has also led 
to the discovery of a much higher degree of complexity 
of GC color-magnitude diagrams. In fact, while some 
clusters seem to photometrically comply with the simple 
stellar population paradigm, a growing number of them 
are found to be characterized by multiple main seque nces 
and/or subgiant and/or giant branches, (jPiotto et al.1 
l2007t IMilone et aLll2008l e.g.,), which have been associ¬ 
ated to variati ons in the content of He and CNO, as we ll 
as age spread (jDAntona et al.ll2005l : iCassisi et al.l l20081. 

This observational evidence led to a general scenario 
where GCs host multiple stellar populations. These are 
often assumed to be associated with different stellar gen¬ 
erations: the ejecta of the slightly older stars, probably 
mixed with varying amounts of gas from the original star 
forming cloud, cr eates subsequent younger generation 
of st ars (see e.g., iDecressin et al.1 [20071 : iD’Ercole et al.1 
2008), though al ternative scenarios are also being con¬ 
sidered (see, e.g.. lBastian et al1 l20131. It is believed that 
only a fraction of the first-generation of stars can con¬ 
tribute to the internal enrichment. The difference in ages 
among the stellar generations is actually relatively small 
for the majority of the clusters (with a few exceptions 
such as e.g., w Cen or M22) and is confined to a couple 
of hundreds Myr. 

The details of this formation scenario are still far from 
being understood. The origin of the polluting material 
remains to be established and it has obvious bearings 
on the timescales for the formation of the cluster it¬ 
self and its mass budget. The observed wide star-to- 
star variations in C, N, O, Na, and A1 found in each 
Galactic GCs, coupled with the uniformity in Fe and Ca 
(apart from a few notable exceptions) provide quite strin¬ 
gent constraints and argue again st anything but a mino r 
contribution from supernovae dCarretta et al.1 l2009alh 
Proposed candidate polluters include intermediate mass 
st ars in t heir asy mptotic giant branch (AGB) phase 
(|Ventura et all [200lh . fast rotating mas sive stars losing 
mass during t heir main sequence phase (IDecressin et al.l 
l200l . novae (iM accaro ne fe Zurekll2012f) and massive bi¬ 
naries (jde Mink et al.ll2009f ). These potential contribu¬ 
tions obviously operate on different time scales and re¬ 
quire a different amount of stellar mass in the first- 
generation. All of the candidates proposed so far fall 
short of reproducing the full variety of observations. Ad¬ 
vances in the theoretical modeling of star formation and 
evolution are likely needed to improve our understanding 
of these issues, including the spanning of a larger range 
of the parameter space (e.g., mass, metallicity, and mass 
loss). However, from an observational point of view, the 
increase in the high-quality abundance work available 
for GCs, both in the sheer number of stars and clus¬ 
ters, as well as in terms of chemical species considered, is 
paramount, as it creates a more complete picture of the 
phenomena involved. 
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The Apache Point Observatory Galactic E volution Ex¬ 
periment (APOGEE iMaiewski et al.l I2015T) is a three 
year, near-infrared (15,090 to 16,990 A; iWilson et al.1 
l2012T >. high-resolution spectroscopic survey of about 
100,000 red giant stars included as part of the 3rd Sloa n 
Digital Sky Survey (SDSS-III lEisenstein et al.l 1201 111 . 
With a nominal resolving power of 22,500, APOGEE 
is deriving abundances of up to 15 elements for nearly 
100,000 stars, although fewer elements are generally de¬ 
tected in weak-lined metal-poor stars. APOGEE is in a 
unique position among the v arious Galactic spect roscopic 
surveys such as Gaia-E SO, dGilmore et al.ll2012f). RA VE 
dSteinmetz et al.l[2006l h and GALAH (lFreemanH2012l) . as 
it uses the_Sloan 2.5m telescope at Apache Point Obser¬ 
vatory (iGunn et alj|2006 ). and thereby has access to the 
northern hemisphere. APOGEE observes a large sample 
of northern globular clusters, something that makes it 
possible to analyse these clusters in a homogeneous way, 
which has not been done before for these objects. 

The study of GCs with APOGEE plays an important 
role not just because it has access to many northern GCs. 
Its high-resolution near-IR spectra allow the simultane¬ 
ous determination of many elemental abundances gen¬ 
erally not available in optical spectroscopic work of GC 
stars. C and N, which are elements heavily affected by 
the pollution phenomenon in GCs, are often not included 
in studies of metal poor stars because the strongest fea¬ 
tures (CH and CN) lie in the near-UV, far from the op¬ 
tical lines of Na, Mg and Al, and thus multiple detectors 
or setups are required to obtain both sets. In addition, 
because these studies usually focus on fairly red stars 
longer exposure times are required to acquire sufficient 
SNR to analyse the near-UV features. 

The spectra used in this paper are publicly ava il able 
as part of the tenth data release (DRIO. lAhn et alll20l4) 
of SDSS-III. The initial set o f star s selected were the 
same used by iMeszaros et al.l (120131 ) to check the accu¬ 
racy and precision of APOGEE parameters published in 
DR10. However, instead of using the automatic ASP- 
CAP pipeline, we will make use of photometry and theo¬ 
retical isochrones to constrain the effective temperature 
(T e g) and surface gravity log g , and use an independent 
semi-automated method for elemental abundance deter¬ 
mination for 10 northern GCs. Some of these clusters are 
well studied, such as M3, M13, M92, M15, while others 
have been poorly studied (NGC 5466), or been only re¬ 
cently discussed in the literature, such as M2 (lYong et al.1 

mi . 

2. ABUNDANCE ANALYSIS 
2.1. Target Selection 

Table 1 lists the globular clusters APOGEE observed 
in its first year, along with the adopted [Fe/H], E(£? —V), 
and ages from the literature. Targets were selected as 
cluster members if 1) there is published abundance infor¬ 
mation on the star as a cluster member, 2) the star is a 
radial velocity member, or 3) if it has a probability >50% 
of being a cluster member based on their proper motion 
adopted from the literature. After this initial selection, 
we checked the position of stars in the T e g — log g dia¬ 
gram based on APOGEE observations and deleted those 
that were not red giant branch (RGB) stars. Stars that 
have metallicity 0.3 dex (typically 3 a scatter) larger or 
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Table 1 

Properties of Clusters from the Literature 


ID 

Name 

N a 

[Fe/H] b 

E(B—V)Ref. c 

NGC 7078 

M15 

23 

-2.37±0.02 

0.10 

NGC 6341 

M92 

47 

-2.31±0.05 

0.02 

NGC 5024 

M53 

16 

-2.10±0.09 

0.02 

NGC 5466 


8 

-I.98i0.09 

0.00 

NGC 6205 

M13 

81 

-1.53±0.04 

0.02 

NGC 7089 

M2 

18 

-1.65±0.07 

0.06 

NGC 5272 

M3 

59 

-1.50±0.05 

0.01 

NGC 5904 

M5 

122 

-1.29±0.02 

0.03 

NGC 6171 

M107 

42 

-1.02±0.02 

0.33 

NGC 6838 

M71 

12 

-0.78±0.02 

0.25 


a N is the number of stars observed in each cluster. 
b [Fe/H] references: IHarris 1996 I <12010 editionl) , clusters are 
listed in order of the average cluster metallicity determined 
in this paper. 

c E(B—V) references: (1) IHarris 1996 I (12010 editionl) . 

smaller than the cluster average also need to be deleted, 
but this last step resulted in no rejections. The clus¬ 
ter target_selection process is described in more detail 
in IZasowski et al.l (l2013l j. The final sample consists of 
428 stars from 10 globular clusters. High S/N spectra 
are essential to determine abundances from atomic and 
molecular features, thus al l selected target s have at least 
S/N=70 as determined bv lMeszaros et all (|2013f ). 

2 . 2 . Atmospheric Parameters 

Abundances presented in this paper are defined for 
each individual element X heavier than helium as 

[X/H] = \og w {n x /n H ) s ^ ~ log w {n x /n H )o (1) 

where n x and uh are respectively the number of atoms 
of element X and hydrogen, per unit volume in the stellar 
photosphere. 

To derive abundances from stellar spectra, we first 
have to estimate four main atmospheric parameters: T e ff, 
log g , microturbulent velocity, and overall metallic¬ 
ity ([Fe/H]). In the following sub-sections we present 
our methodology for determining these parameters and 
our reasons for not using the values available for each 
star from the APOGEE Stellar Para meters and C hemi¬ 
cal Abundances Pipeline (ASPCAP; iGarcia Perez et all 
l2014fl in DR10. In order to evaluate th e accuracy and 
precis ion of the ASPCAP parameters, Mcszaro s et all 
(12013H carried out a careful comparison with literature 
values using 559 stars in 20 open and globular clusters. 

These clusters were chosen to cover most of the pa- 
rameter ra nge o f star s APOGEE is expected to observe. 
iMeszaros et al.l (|2013l f provided a detailed explanation of 
the accuracy and precision of these parameters, and also 
derived empirical calibrations for T e ff, log g , and [M/H] 
using literature data. In the sections below we will briefly 
review these calibrations along with their limitations. 

2 . 3 . The Effective Temperature 

We adopt a photometric effective temperatures cal- 
cula ted from the J — K s colors using the equations of 
IGonzalez Hernandez fe Bonifaciol (|2009l l. Their calibra¬ 
tion was chosen because of its proximity of only 30—40 K 
to the absolute temperature scale. De-reddened J — K s 


were calculated the same way as by [Mcszaro s et al.1 
(l2013i l. from E(B — V ), listed in Table 1 for each cluster, 
using E(J - K s ) = 0.46 • E(B - V). 

The ASPCAP DR10 effective temperatures were 
com pared with photo metr ic o nes u sing _calibrations 
by IGonzalez Hernandez fe Bonifaciol (120090 base d 
on 2MASS J — K„ colors dStrutskie et all 120061) . 
IMeszaros et al.l (120131) found that small systematic 
differences, in the range of 100—200K, are present 
between the ASPCAP and photometric tempera¬ 
tures. The ASPCAP DR10 T e g were also compared 
to literature spectroscopic temperatures, and the 
average of these differences were found to be negli¬ 
gible. The corrected ASPCAP DR10 temperatures 
were calculated between 3500 K and 5500 K using 

a calibration relati on derived from_the comparison 

with the IGonzalez Hernandez fe Bonifaciol (1200911 scale. 
ASPCAP DR10 raw temperatures above 5000 K in 
metal-poor stars showed significant, 300—500 K offsets 
compared to photometry, and are thus believed to be 
not accurate enough for abundance analysis. This issue 
is mostly limited to metal poor stars, and does not 
affect stars at [Fe/H]> —1, where the vast majority of 
APOGEE targets are. 

The adoption of a purely photometric temperature 
scale enables us to be somewhat independent of ASP¬ 
CAP (while still using the same spectra), which gives im¬ 
portant comparison data for future ASPCAP validation. 
Besides providing an independent comparison dataset for 
APOGEE, the photometric temperatures allowed us to 
include stars in the sample that are hotter than 5000 K. 
Because of these reasons, the final results presented in 
this paper are based on the photometric temperatures, 
and we only use the ASPCAP DR10 raw temperatures 
to estimate our errors related to the atmospheric param¬ 
eters. 

2 . 4 . The Metallicity 

The APOGEE DR10 release contains metallicities de¬ 
rived by ASPCAP for all stars analysed, thus provid¬ 
ing an alternative scale to manually derived metallici¬ 
ties. That metallicity, [M/H], tracks all metals relative 
to the Sun, and gives the overall metallicity of the stars 
because it was derived by fitting the entire wavelength 
region covered by the APOGEE spectrograph. This is 
different from most literature publications that use Fe 
lines to track metallicity in a stellar atmosphere. We use 
[Fe/H] in this paper whenever we refer to metallicity pre¬ 
sented here, because we use Fe I lines to measure it. For 
the most part, one can treat values of [M/H] as if they 
were [Fe/H]. When ASPCAP metallicity was compared 
with individual values from high-resolution observations 
from the literature, a difference of 0.1 dex is found below 
[M/H]=—1, and this discrepancy increased with decreas- 
ing metallicity reachin g 0.2—0.3 dex around [M/H] =—2 
(IMeszaros et all [201311 . The calibrated DR10 metallici¬ 
ties map well onto [Fe/H], because the calibration process 
uses [Fe/H] values from the optical. 

Because an alternative metallicity based on only iron 
lines did not exist for APOGEE in DR10, and because 
of these small zero point offsets, we decided to derive our 
own [Fe/H] based on Fe I lines found in the H band, in¬ 
stead of using the published APOGEE DR10 [M/H] val¬ 
ues. We assumed that all stars have the literature clus- 

























































4 


Table 2 

Properties of Stars Analyzed 


2MASS ID 

Cluster 

v helio 

T e fl 

log g 

[Fe/H] 

[C/Fe] 

[N/Fe] 

[O/Fe] 

[Mg/Fe] 

[Al/Fe] 

[Si/Fe] 

[Ca/Fe] [Ti/Fe] 

2M21301565+1208229 

M 15 

-104.5 

4836 

1.56 

-2.12 




0.16 

-0.06 

0.35 

0.19 

2M21301606+1213342 

mis 

-108.3 

4870 

1.64 

-2.31 




0.10 

0.57 

0.46 

0.53 

2M21304412+1211226 

M 15 

-102.7 

4715 

1.28 

-2.12 




-0.45 

0.63 

0.60 

0.35 

2M21290843+1209118 

M 15 

-106.0 

4607 

1.03 

-2.07 




-0.11 

0.75 

0.41 


2M21294979+1211058 

mis 

-107.6 

4375 

0.56 

-2.31 

-0.44 

0.95 

0.44 

0.17 

0.64 

0.44 

0.06 


Note. — This table is available in its entirety in machine-readable form in the online journal. A portion is shown here for guidance regarding its form 
and content. 


ter mean metallicity before starting the calculations, af¬ 
ter which small wavelength windows around the Fe lines 
listed in Table 3 were used to revise the individual star 
metallicities. 


2.5. The Surface Gravity 

In this study we ado pt gravities from stella r evolution 
calculations. Following iMeszaros et al. ( 201 3b we derive 
graviti es for our sample using isoc hrones from the Padova 
group (|Bertelli et all I2008L120091) . The cluster metallic¬ 
ities collected from the literature used in the isochrones 
are listed in Table 1, while the ages of all clusters were 
chosen to be 10 Gyr. The final set of temperatures and 
gravities corresponding to them from the isochrones are 
listed in Table 2. 

The ASPCAP DR10 surface gravit ies were com- 
pared to both isochrones and Kepler dBorucki et al. I 
l2010f) asteroseismic targ ets observed by APOGEE 
( Pinsonneault. et al.1120141 ) to estimate their accuracy 
and precision. An average difference of 0.3 dex was 
found at solar metallicity in both cases, but this in¬ 
creased to almost 1 dex for very metal poor stars. The 
asteroseismic gravities are believed to have errors in the 
range of 0.01—0.05 dex, thus far superior to spectro¬ 
scopic measurements, but they are only available for 
metal rich stars with [M/H]> —1.0. The discrepancy be¬ 
tween spectroscopic and asteroseism ic su r face gravities is 
a top ic of ongoing in vesti gation (e.g JPinsonneault. et al.1 
2014 : [Epstein et all 120141 ). The final calibration from 


Meszaros et al.l ( 20131) combined surface gravities derived 


from isochrones below [M/H]< —1.0 with the asteroseis¬ 
mic dataset for high metallicities. Values after the cali¬ 
bration still show some small (<0.1 dex) offsets and large 
scatter around the isochrones for the globular clusters. 
Thus, in order to minimize errors in abundances related 
to the uncertainty in the ASPCAP corrected spectro¬ 
scopic gravities, we decided to use the pure isochrone 
gravities in this study. 


2.6. The microturbulent velocity 

A relation between microturbulent velocity and surface 
gravity, v m i cro = 2 24 — 0.3 x log q, was u sed for ASPCAP 
analyses in DR10 (IMeszaros et al.l[2013l ). For consistency 
with ASPCAP, and for easier future comparisons with 
APOGEE results, we adopted that equation in this work. 

3. AUTO SYNTH 

The program called autosynth was developed especially 
for this project to simplify the large amount of synthe¬ 
sis required for abundance determination. The program 
takes atmospheric parameters as input and carries out 


spectral synthesis to derive elemental _abundances. The 
core of autosynth is MOOG2013 dSnedenl I1973D . which 
does the spectrum synthesis, while autosynth compares 
the synthetic spectrum with the observed spectrum and 
determines the best abundances with y 2 minimization in 
wavelength windows specified by the user. The program 
can read the MOOG formatted ATLAS and MARCS 
model atmospheres, and can also convert the original 
ATLAS and MARCS formats into a MOOG compatible 
format. 

The line list adopted for this study includes both 
atomic and molecular species. It is an updated ver¬ 
sion of what was used for the DR10 results, version 
m201312160900 (als o used for DR12 APOGEE’s next 
publi c data release; lAlam et, al.l 120151 iHoltzman et al.1 
2015), and includes atomic and molecular species. The 
molecular line list is a compilation of literature sources 
including transitions of CO, OH, CN, C 2 , H 2 , and SiH. 
All molecular data are adopted without change with the 
exception of a few obvious typographical corrections. 
The atomic line list was compiled from a number of liter¬ 
ature sources and includes theoretical, astrophysical and 
laboratory oscillator strength values. These literature 
line positions, oscillator strengths, and damping values 
were allowed to vary in order to fit to the solar spectrum 
and the spectrum of Arcturus, thus generating a tuned 
astrophysical line list. The solution is weighted such that 
the solar solution has twice the weight as the Arcturus 
solution to properly consider the fact that the abundance 
ratios in Arcturus are more poorly understood than those 
of the Sun. The code used for this process was ba sed on 
the LTE spectral synthesis code MOOG dSnedenl Il973i) 
but adapted to our unique needs. For lines with labora¬ 
tory oscillator strengths, we did not allow the astrophys¬ 
ical gf value to vary beyond twice the error quoted by 
the source. A more detailed de scripti on of this process 
and the line list can be found in ISlretrone et al.1 (|2015l ) . 

The choice for the local continuum set can greatly af¬ 
fect the derived abundances, thus we needed a reliable 
automated way to determine the continuum placement. 
This was done with a separate y 2 minimization from 
the one that was used for the abundance determina¬ 
tion. Continuum normalized observation points around 
1 are multiplied by a factor between 0.7 and 1.1 with 
0.001 steps for each synthesis emulating slightly different 
choices for the location of the local continuum. Multipli¬ 
cation is necessary because it preserves the original spec¬ 
trum. The x near the continuum is calculated and com¬ 
pared to the continuum of the observation and minimized 

21 http://www.as.utexas.edu/ chris/moog.html 
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separately for every abundance step. The y 2 calculation 
for the abundances determination happens between cer¬ 
tain flux ranges (usually between 0.3 and 1.1) using the 
continuum placement determined in the previous step. 

3 . 1 . Individual Abundances 

The individual abundances were determined using the 
ID Local Thermodynamic Equilibrium (LTE) model at¬ 
mospheres calculated with ATLAS9 (lKuruc^ll979 ). The 
model atmospher es were generated us ing solar reference 
abundances from lAsnlund et ak i (120051 ). the same way as 
the main AP OGEE m o del at mosphere database was gen¬ 
erated dMeszaros et al.ll2012lh Because the overall metal- 
licity of these clusters were well known from the litera¬ 
ture, initially we calculated atmospheres using the av¬ 
erage literature metallicity for each cluster adopting the 
photometric effective temperatures and isochrone gravi¬ 
ties. These initial model atmospheres were later revised 
to have consistency with the synthesis. 

The windows used to derive the individual abundances 
were determined based on the analysis of F TS stars i n 
t he H —band using the APOGEE line list bv lSmitlr et al.l 
( 2013 1. In the case of Fe we measured [Fe/M], relative to 
the literature cluster metallicity for each line. The abun¬ 
dance of Na is very important in discussing the spread of 
O in GCs, and two Na lines are available in the APOGEE 
spectral band. However, these two Na lines are weak even 
at solar metallicities. We carried out a number of tests 
attempting to derive Na abundances, but we found that 
the two Na lines become very weak around [Fe/H]=—0.5, 
and non-detectable below about —0.7, thus we were not 
able to determine Na abundances for any of our targets. 
The list of wavelength regions used in our analysis, and 
the solar reference values for each element is listed in Ta¬ 
ble 3. Figures [T] and [5] show examples of observed Fe, Mg, 
Al, OH, CO, and CN line profiles and their fitted syn¬ 
thesis for one star from M71 and M13. The wavelength 
regions shown in these figures are only a fraction of what 
has been used from Table 3. 

CN lines spread over most of the H—band, hence it is 
important to calculate the CNO abundances before the 
atomic ones. It is also important to use self consistent 
model atmospheres, because stars in globular clusters ex¬ 
hibit low carbon and high a—content, which significantly 
alters the struct ure of the atmosphere compared to a so¬ 
lar scaled one (Meszaros et al. 2012). Taking into ac¬ 
count all this, we developed the following procedure to 
produce the final abundances for each star: 

1. A model atmosphere is generated using literature 
cluster average metallicities, the photometric tem¬ 
perature and an isochrone gravity. Because all 
of our targets are RGB stars, we choose [C/Fe] 
= —0.5, [O/Fe] = 0.3, and [N/Fe] = 0.5 dex for 
this initial model. 

2. Individual Fe I lines are fit with autosynth , and an 
average [Fe/H] is calculated for each star. 

3. A new model atmosphere is calculated using this 
new [Fe/H] value, but still using the starting CNO 
abundances. 

4. We set the abundances of C, N and O before the 
remaining elements, because they can have a sig- 


Table 3 

Wavelength Regions 


Element 

l°g(N) a 

Wavelength (A) b 

Fe 

7.45 

15210-15213.5; 15397-15401; 15651-15654 



15966-15973; 

16044-16048; 16156-16160 
16168-16171 

C 

8.39 

15572-15606; 

15772-15791; 15980-16037 



16172-16248; 

16617-16677; 16839-16870 

N 

7.78 


15240-15417 

O 

8.66 

15267-15272; 

15281-15288; 15372-15380 



15386-15390; 

15394-15397; 15404-15414 



15499-15502; 

15508-15511; 15539-15542 



15561-15566; 

15569-15574; 15887-15904 



16188-16198; 

16207-16213; 16233-16237 



16244-16247; 

16251-16261; 16300-16305 



16314-16319; 

16707-16714; 16718-16720 



16731-16735; 

16888-16892; 16898-16912 

Mg 

7.53 

15741- 

15757; 15767-15773 

Al 

6.37 

16720-16727; 

16751-16759; 16765-16770 

Si 

7.51 

15962-15966; 

16062-16066; 16097-16101 



16218-16223; 

16683-16687; 16830-16834 

Ca 

6.31 

16139-16143; 16153.5-16164 

Ti 

4.90 

15546.5-15549.5; 15718-15721.5 


a The Solar reference abundances are from lAsplund et al] r200."ii l. 
b Vacuum wavelength. 


nificant effect on the atmospheric structure in cool 
stars. Since molecular features generally disappear 
from metal-poor spectra above 4500 K, we divide 
our stars into two temperature groups. For the 
stars cooler than 4500 K, we first determine [O/Fe] 
using OH lines, then create a new model atmo¬ 
sphere with [a/Fe] equal to [O/Fe]. We then de¬ 
termine C and O abundances from CO lines, then 
recreate the model atmosphere again with these 
new [C/Fe] and [O/Fe] abundances. Finally, we de¬ 
rive N abundance using CN lines. For stars hotter 
than 4500 K, we leave the C, N, and O abundances 
at their inital values. 

5. The abundances of the remaining elements (Mg, 
Al, Si, Ca and Ti) are determined with autosynth, 
using the stellar parameters, metallicities, and C, 
N and O abundances previously determined. 

For each element, we average together the abundance 
results from the different wavelength regions to obtain 
final values. Although the size of each region is differ¬ 
ent, we did not find it necessary to use weights based 
on their ranges or line strengths, because that approach 
did not produce abundances significantly different from a 
straightforward average. Data reduction errors or miss¬ 
ing data affected some of these regions, resulting in er¬ 
roneous fits, and because of this we carefully examined 
each fit by eye. These wavelength regions were not in¬ 
cluded when constructing the final average abundances. 
The final abundance values are listed in Table 2. 

3 . 2 . UnceHainty Calculations 
3 . 2 . 1 . Systematic Uncertainties 

The uncertainty in the atmospheric parameters 
strongly affects the final abundances derived from some 
of the spectral features we consider. To test the sensi¬ 
tivity of abundances due to changes in the atmospheric 
parameters we used the results from the ASPCAP raw 
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Figure 1 . Example spectra and the fitted synthesis of Fe, Al, and Mg lines for two stars from M71 and M13. Abundances were fitted 
between the labelled minimum and maximum values using a step of 0.01 dex. The printed best fitted abundance values might not be the 
same as in Table 2, because the table contains averaged values, not individual fits. 


temperature scale. The same exact steps described in 
the previous section were followed, but instead of adopt¬ 
ing the photometric temperature scale we adopt the AS- 
PCAP DR10 raw temperature scale, which results in 
new surface gravities and microturbulent velocities. This 
way we could track systematics uncertainties sensitive to 
these parameters as well. 

The differences in abundances as a function of photo¬ 
metric temperatures are demonstrated in Figure [3] The 
top left panel displays the differences in the measured 
abundances by using ASPCAP and photometric temper¬ 


ature, while the rest of the panels are assigned to each 
element. The color scale in all panels represents A T e ff. 
We defined the estimated errors associated with the at¬ 
mospheric parameters based on the standard deviation 
around the mean differences between the two tempera¬ 
ture scales. The calculated standard deviation in of the 
differences in temperatures is 146 K (which we round to 
150 K). This standard deviation corresponds to the sum 
of the uncertainty in the photometric temperature and 
the ASPCAP temperature in quadrature. 

In order to estimate the uncertainty of just the photo- 
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L 
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Wavelength (A) 


Figure 2. Example spectra and the fitted synthesis of OH, CO, and CN lines for two stars from M71 and M13. For more explanation see 
caption of Figure [l] and Section 3.1. 


metric temperature component, we carried out calcula¬ 
tions of temperatures with varied J—Ks colors, redden¬ 
ings and metallicities for M107. M107 was chosen be¬ 
cause it is the cluster with the highest reddening in our 
sample. The uncertainty of 2MASS photometry is usu¬ 
ally between 0.025—0.03 magnitude for these stars, and 
by using 0.03 as baseline, we estimate an uncertainty of 
0.05 magnitude for the J—Ks color. Changing J—Ks by 
0.05 magnitudes typically produces a change of 80 K in 
the photometric temperature. The reddening of M107 
was changed by 0.03, simulating a 10% uncertainty in 
reddening, and this produced a difference of about 40 K 


in the photometric temperature. A change of 0.1 dex 
in metallicity results in 1 K or less uncertainty in tem¬ 
perature; thus the uncertainties in metallicity can be ne¬ 
glected. By adding the uncertainty from photometry and 
reddening in quadrature, we estimate the uncertainty of 
the photometric temperature to be about 100 K. The top 
panel of Table 4 lists uncertainties associated with both 
150 and 100 K changes in temperature. 

In our methodology a 100 K change in temperature 
also introduces an 0.3 dex systematic difference in sur¬ 
face gravity, and 0.1 km/s in microturbulent velocity, so 
by coupling these two parameters to T e ff, our system- 



































































Table 4 

Estimated Abundance Uncertainties 


Cluster Fe C N O Mg A1 Si Ca Ti 


Systematic uncertainties from atmospheric parameters 


AT ef f 
= 150 I < 

0.08 

0.12 

0.16 

0.18 

0.04 

0.07 

0.05 

0.06 

0.15 

ATgff 
= 100 I < 

0.05 

0.08 

0.11 

0.12 

0.03 

0.05 

0.03 

0.04 

0.10 


Random uncertainties from autosynth 


M 15 

0.11 

0.21 

0.30 

0.06 

0.08 

0.16 

0.13 

0.19 

0.26 

M 92 

0.10 

0.18 

0.30 

0.05 

0.07 

0.12 

0.12 

0.22 

0.17 

M 53 

0.06 

0.16 

0.12 

0.05 

0.06 

0.14 

0.11 

0.11 

0.18 

N 5466 

0.06 

0.18 

0.15 

0.03 

0.07 

0.03 

0.07 

0.14 

0.17 

M 13 

0.05 

0.12 

0.06 

0.05 

0.05 

0.12 

0.07 

0.06 

0.12 

M 2 

0.04 

0.11 

0.06 

0.04 

0.07 

0.10 

0.09 

0.05 

0.09 

M 3 

0.04 

0.12 

0.08 

0.05 

0.03 

0.10 

0.08 

0.05 

0.12 

M 5 

0.04 

0.08 

0.05 

0.05 

0.04 

0.06 

0.07 

0.04 

0.09 

M 107 

0.04 

0.04 

0.03 

0.05 

0.04 

0.08 

0.07 

0.04 

0.10 

M 71 

0.02 

0.03 

0.04 

0.04 

0.04 

0.05 

0.07 

0.06 

0.07 


Final combined uncertainties 


M 15 

0.12 

0.22 

0.32 

0.13 

0.09 

0.17 

0.13 

0.19 

0.28 

M 92 

0.11 

0.20 

0.32 

0.13 

0.09 

0.13 

0.12 

0.22 

0.20 

M 53 

0.08 

0.18 

0.16 

0.13 

0.07 

0.15 

0.11 

0.12 

0.21 

N 5466 

0.08 

0.20 

0.19 

0.12 

0.08 

0.06 

0.08 

0.15 

0.20 

M 13 

0.07 

0.14 

0.13 

0.13 

0.06 

0.13 

0.08 

0.07 

0.16 

M 2 

0.06 

0.14 

0.13 

0.13 

0.08 

0.11 

0.09 

0.06 

0.13 

M 3 

0.06 

0.14 

0.14 

0.13 

0.04 

0.11 

0.09 

0.06 

0.16 

M 5 

0.06 

0.11 

0.12 

0.13 

0.05 

0.08 

0.08 

0.06 

0.13 

M 107 

0.06 

0.09 

0.11 

0.13 

0.05 

0.09 

0.08 

0.06 

0.14 

M 71 

0.05 

0.09 

0.12 

0.13 

0.05 

0.07 

0.08 

0.07 

0.12 


Note. — Top panel: Systematic uncertainty estimates from changes in 
T e ff, log g, and u m i cro . Middle panel: Average random uncertainties re¬ 
ported by autosynth. Bottom panel: The final uncertainties are the sum 
of uncertainties in quadrature from the middle panel and uncertainties for 
AT e ff = 100K, A log p=0.3, and A Vmicro— 0.1 km/s. 

atic uncertainties include the investigation of abundance 
sensitivity to these parameters as well. 

3.2.2. Internal Uncertainties 
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Figure 3. Differences in abundances produced by two runs adopt¬ 
ing different temperatures: photometric and ASPCAP tempera¬ 
tures; otherwise the same calculation method was used. The points 
are color coded by the differences between the photometric and 
ASPCAP temperatures. The ± errors give the standard deviation 
around the mean of the differences. 


Besides the uncertainty coming from the adopted at¬ 
mospheric parameters, the uncertainty of the fit was also 
calculated using the a of the residuals between the best 
fit synthesis and the observation. These calculations esti¬ 
mate random uncertainties. This a of the fit is calculated 
within the windows used in the y 2 calculation. For de¬ 
termining the uncertainty of the fit, we multiply the ob¬ 
served spectra by 1+cr and l-cr which simulates two spec¬ 
tra slightly different from the original spectrum. Then, 
the fit of each line is repeated while keeping all other 
parameters unchanged. The average differences between 
these two new fits and the original best-fit spectrum is 
the defined uncertainty associated with the fit itself. This 
uncertainty estimate is mainly sensitive to variations of 
noise in the spectrum in the defined windows for the x 
fit. If its value is close to 0, then the (1 + cr)x spectrum 
and (1 — a) x spectrum give very similar, or the same 
abundances. This is expected when working at high S/N 
and with a well defined continuum. 

While this uncertainty estimation method is reliable 
in most cases, it has its limitations for very noisy spec¬ 
tra, very weak lines, or when abundances are near upper 
limits. Uncertainties are usually overestimated for noisy 
spectra, while they are underestimated for weak lines 


and upper limits. Thus, we decided to use one uncer¬ 
tainty from autosynth per element per cluster by simply 
averaging together all uncertainties reported by the pro¬ 
gram. This resulted in one uncertainty estimate for ev¬ 
ery element in each cluster. Autosynth uncertainties are 
listed in the middle section of Table 4. The final esti¬ 
mated uncertainties were calculated by adding together 
in quadrature the uncertainties for 100 K difference in 
temperature (also 0.3 dex systematic difference in log g 
, and 0.1 km/s in microturbulent velocity), and the aver¬ 
age autosynth values that estimate random uncertainties. 
These final estimations for each element per cluster are 
given in the bottom panel of Table 4. 

4. LITERATURE COMPARISONS 

We take [X/Fe] abundance values for C, N, O, Mg, Al, 
Si, Ca and Ti from high-resolution spectroscopic stud¬ 
ies in the literature as a point of comparison for our 
abundance determinations . We u se the same literature 
sources as iMeszaros et al. (201 3). and added more re¬ 
cently published papers listed in Table 5. The deriva¬ 
tion of stellar parameters T e ff and log g are described 
in detail in Section 2; and they were also compared to 
the literature. Our goal for these comparisons is not to 
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Table 5 

Identifiers, stellar parameters and elemental abundances from the literature for stars in our sample 


2MASS ID Cluster T eff log g [Fe/H] [C/Fe] [N/Fe] 


2M21290843+1209118 M15 -2.37 

2M21293353+1204552 M15 -2.37 

2M21293871+1211530 M15 -2.37 

2M21294979+1211058 M15 -2.40 

2M21294979+1211058 M15 -2.27 


[O/Fe] [Mg/Fe] [Al/Fe] [Si/Fe] [Ca/Fe] [Ti/Fe] alt. ID Sour 


0.42 


0.36 


0.31 B5 

h 


0.13 

K22 

h 


0.14 

K47 

h 

0.53 

0.28 

K144 

h 


0.38 

K144 

i 


Note. 


This table is avai lable in its entire t y in m achi ne-readable form in the onl ine journal. A portion is shown her e for guidance regarding i ts form and contei 

(c)l Joh nson_fc_Pilachowski ( 2012ibldl|Sneden etj+ N2004l :/+|Yonj! L etaU)|2006; 

(llj^^nrh^e^^"Tlfl96[k rnjOtsi]^r^t*an~il200l|h n£lJpnedeimj^^J^i99l|h 


Individual star referenc es: (a) [CVComie^^^^^J2^HJ) ; (b) Carrctta et al.| l|2009a| l: fel l 
(f) Cohen V Melendez (2005); (g) Cavallo V Nagar (2000); (h) Sncden et ah, (2000) 


I 5:ii( :* •: w*T77"]~]TtT^ 7 1 nil |S()i)e('k e ~'"TH~TbG l I l| : '; ]KT77t‘v,Iv;mT /20()((|: 1 : ((O' [Kra(T~/t aI.| :( ^92 


{ 2 Oil]) ; (r) 


1: f s i lKoc^^^Mc ^/ilUamn 201t]jh_itT[Sn^^n_et al.l TTilO^lT Ini I ilainire/ fc~CohR^j200^h_fvI Yong_^^TT l200^ b_lw ) \ Ic le 11 dt TT^^Joh en] f200T: (x) mril^o 
b_TvVTsh^ro^^^9^i: Izl ISmith et g^ ~ ii200^1 : (aa ULee et al] joQQil) : !abT ^b^I^ 2 ^z Cohen i 200(3): lari lYona et g^ " iT2006bl' ) : ladT TSneden et, ah! 120001 : 
Roederer & Sncden (2011) 


cross-calibrate our new abundance determinations with 
the literature; rather, we are looking for cases where our 
abundances are systematically different from the litera¬ 
ture, or particular clusters or elements for which our ho¬ 
mogeneously observed and analysed data set can clarify 
conflicts in the literature. 

Cross-identification between the globular cluster stars 
in the DR10 APOGEE data set and the literature was 
performed using the Simbad online serviciF^I. based on 
2MASS identifiers and (a, 6) coordinates. Because there 
is a large and heterogeneous literature on chemical abun¬ 
dances in globular cluster stars, we are providing our 
cross-identifications as a resource for the community. Ta¬ 
ble 5 lists 2MASS ID and position, and alternate stellar 
identifiers from literature abundance studies for all of the 
stars considered in this study. 

Figure []] shows comparisons of our stellar parameters 
and abundances against the literature values. Differ¬ 
ent globular clusters are represented by different colored 
points. In general, we find a systematic offset of ~ 100 
K between our photometric T e ff values and the spectro¬ 
scopic effective temperatures from the literature, with a 
reasonably small scatter. This is similar to the typical 
difference between spectr osco p ic and photometric tem¬ 
peratures reported by IMeszaros et all (120131) . Because 
of the degeneracies in deriving stellar parameters, the 
slightly higher temperatures in the literature are accom¬ 
panied by a systematic offset of ~ +0.2 dex in log g. 

There are a few systematic differences between our 
abundance results and those in the literature. These can 
mainly be traced back to a change in the S olar abundance 
scale as derived bv lAsplund et al.l (120051) . As one exam¬ 
ple, the [Fe/H] metallicities we derive are typically 0.1 
dex higher than those from the literature, which is quite 
similar in magnitude to th e chang e in A(Fe)^ from 7.5 2 
(|Anders &; G re vessel 1 989Tl to 7.45 (lAsplund et al.ll2005l) . 
Also, the Solar abundance of oxygen was revised s ignif¬ 
icantl y, from A(Q)f7i=8.93 (lAnders &: G revessel 1 19891) to 
8.66 (lAsplund et a. 1.1 [20051) . Since we use the m ore re - 
cent Solar abundance values from lAsplund et al.l (120051 ) 
whereas our earlier literature sources do not, we expect 
the ~ 0.3 dex offset visible in Figure 0] The systematic 
difference in [N/Fe] abundance is also likely to be due to 
the updated Solar abundances. 

Carbon is the only element studied in this paper for 

22 http: //simbad.u-strasbg.fr/simbad/ 


which this explanation does not hold: although the 
Solar abundance was rev ised dow n, from A(CW=8.5 6 
(lAnders fe Grevesselll989D to 8.39 (lAsplund et alJl2005h . 
our [C/Fe] values are typically lower than the literature 
values. Unfortunately, literature values are available for 
only a subset of our stars, which makes it difficult to 
verify precisely any systematic behavior. 

In the other elements we consider in the present study, 
there is generally good agreement between our results 
and those from the literature, which is encouraging. 
However, we find significant differences in Al abundances 
in stars with below-Solar metallicities. This is mainly 
driven by stars i n M3, where the [Al/Fe] abundances 
from iJohnson et al.l (120051) are larger than in the other 
literature sources. We also see significant scatter in 
[Ca/Fe] and [Ti/Fe] at low abundances; this is likely 
caused by our abundance determination method having 
difficulty with weak lines. 

The [Mg/Fe] versus [Al/Fe] relations in M3, M13 and 
M5 provide useful examples of how our new data set 
compares with the literature. In M3, our abundance de¬ 
terminations show a clea r bimodality in the Mg- Al dis¬ 
tribution. The study of ICavallo fe Nagail (120001) found 
a similar bimodality, but IJohnson et al.l (120051) found a 
smooth distribution. In M13 we find an extended an¬ 
ticorrelation between [Mg/Fe] and [Al/Fe], and the lit¬ 
erature abunda nces h e within it. We have seven stars 
in common with iSneden et al.l (120041 ), and they span the 
full anticorrelation, while the three a nd two stars, respec¬ 
tively, that we have in comm o n withICohen fe Melendez! 
(120051) and ICavallo fe Nagail (|2000l ) are consistent with 
our abundance results but happen to inhabit small re¬ 
gions of abundance space. M5 is a similar case, where 
our [Al/Fe] val ues span a range of —0.3 dex to +1.1 dex, 
while the ICarretta et al.l (12009a ) Al abundances are lim¬ 
ited to between —0.2 and +0.7 dex. As in M13, our 
sample is larger than the literature sample, and covers 
the full RGB, while the literature sample does not fully 
span the parameter space. 

5. VARIATIONS IN INDIVIDUAL ELEMENT ABUNDANCES 

Although there are well-known abundance patterns 
within globular clusters, large homogeneous studies that 
include a wide range of abundances are rare. Since we 
can determine abundances of most of the light elements 
for the stars in our sample, we examine the behavior 
of the known abundance patterns across a wide range 
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Literature Comparisons 



5000 4500 4000 


2.5 2 1.5 1 0.5 0 
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[C/Fel this study 






Figure 4. T e ff, log g, [Fe/H], [C/Fe], [N/Fe], [Mg/Fe], [Al/Fe], [Si/Fe], [Ca/Fe], [Ti/Fe] compared to literature sources. Different colors 
denote different clusters, for an explanation of the colors see upper right panel. Solid lines show the 1:1 relation, while dashed lines denote 
±100 K for T e ff, and ±0.2 dex for log g and individual abundances. A detailed discussion can be found in Section 4. 


in cluster metallicity, and search for unexpected varia¬ 
tions in the a—elements. In this section we focus on the 
range of abundance within each cluster, to separate real 
abundance variations, bimodalities and trends from pos¬ 
sible measurement errors. Table 6 lists the average and 
standard deviation for each elemental abundance in each 
cluster. 

5.1. Correlations with T e g 

In Figures 07 ] we show all nine derived abundances 
as a function of effective temperature. From Figure [5] 
we conclude that we measure constant Fe abundances in 


all of the clusters. Mg and A1 abundances show a large 
range of values in some clusters (as discussed in Section 
2, and further in Section 6), but no significant trends 
with T e ff except in M107 and M71. This trend is very 
weak in M71, and data are consistent with showing no 
trend within the uncertainties. However, in M107 the 
trend is stronger. We currently do not fully understand 
where these small correlations come from. We suspect 
that this is a result of a combination of effects. One 
such effect may be the use of model atmospheres that as¬ 
sume local thermodynamic equilibrium (LTE), but that 
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Table 6 

Abundance Averages and Scatter 


Cluster 

[Fe/H] 

[C/Fe] 

[N/Fe] 

[O/Fe] 

A(C+N+0) 

[Mg/Fe] 

[Al/Fe] 

[Si/Fe] 

[Ca/Fe] 

[Ti/Fe] 


Averages 

M15 

-2.28 

-0.41 

0.95 

0.54 

7.09 

0.11 

0.34 

0.44 

0.16 

0.15 

M92 

-2.23 

-0.41 

0.93 

0.58 

7.19 

0.14 

0.42 

0.45 

0.10 

0.09 

M53 

-1.95 

-0.50 

1.06 

0.56 

7.49 

0.11 

0.37 

0.41 

0.23 

0.28 

N5466 

-1.82 

-0.56 

0.84 

0.63 

7.60 

0.14 

-0.24 

0.29 

0.04 

0.29 

M13 

-1.50 

-0.53 

0.89 

0.28 

7.69 

0.13 

0.61 

0.40 

0.26 

0.20 

M2 

-1.49 

-0.48 

0.90 

0.41 

7.76 

0.26 

0.45 

0.35 

0.24 

0.27 

M3 

-1.40 

-0.46 

0.69 

0.40 

7.84 

0.15 

0.21 

0.30 

0.12 

0.11 

M5 

-1.24 

-0.46 

0.76 

0.27 

7.85 

0.23 

0.36 

0.34 

0.20 

0.26 

M107 

-1.01 

-0.21 

0.69 

0.33 

8.15 

0.24 

0.47 

0.48 

0.15 

0.21 

M71 

-0.68 

-0.10 

0.91 

0.51 

8.65 

0.38 

0.51 

0.39 

0.21 

0.42 

Scatter 

M15 

0.10 

0.13 

0.35 

0.19 

0.14 

0.24 

0.52 

0.16 

0.25 

0.08 

M92 

0.10 

0.11 

0.23 

0.19 

0.13 

0.23 

0.48 

0.12 

0.17 

0.20 

M53 

0.07 

0.16 

0.21 

0.06 

0.09 

0.08 

0.51 

0.05 

0.17 

0.13 

N5466 

0.08 

0.01 

0.10 

0.04 

0.08 

0.06 

0.35 

0.09 

0.25 

0.14 

M13 

0.07 

0.07 a 

0.17 

0.17 

0.15 

0.15 

0.53 

0.09 

0.14 

0.14 

M2 

0.08 

0.05 a 

0.15 

0.21 

0.23 

0.07 

0.51 

0.06 

0.12 

0.08 a 

M3 

0.08 

0.08 a 

0.13 

0.23 

0.15 

0.06 

0.43 

0.04 

0.08 

0.13 

M5 

0.08 

0.13 a 

0.17 

0.27 

0.19 

0.09 

0.34 

0.07 

0.11 

0 .12 a 

M107 

0.06 

0.09 

0.27 

0.15 

0.15 

0 .10 a 

0.15 

0.05 

0.21 

0.16 

M71 

0.07 

0.13 

0.32 

0.09 

0.13 

0.05 a 

0.08 

0.05 

0.10 

0.09 


a Standard deviation around the linear fit. 


[Fe/H], [Mg/Fe], [Al/Fe] as a function of T eff , T eff < 5400K 
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Figure 5. [Fe/H], [Mg/Fe], and [Al/Fe] as a function of photometric T e ff for all ten clusters. The error bars represent our final combined 
uncertainties from Table 4. The linear fit is plotted over [Mg/Fe] for M107 and M71 to remove the visible linear trend (see Section 5.1 for 
discussion). 

non-LTE effects (iBergemann fe Nordlanderll2014f) act to make the strong Mg lines in these metal-rich stars give 
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Figure 6. [C/Fe], [N/Fe], and [O/Fe] as a function of photometric T e fj for all ten clusters. Open triangles mark upper limits for [C/Fe] and 
lower limits for [N/Fe], while the real detections are plotted using filled red dots. The error bars represent our final combined uncertainties 
from Table 4. The linear correlation in [C/Fe] as a function of T e ff in M13, M2, M3 and M5 is the effect of CNO burning on the RGB. 
The fitted lines are used to remove the trend in order to estimate the scatter in these clusters (see Section 5.1 for discussion). 


the appearance of higher abundance in the cooler stars. 
Other possible effects are small systematic errors from 
estimating T e g and log g culminating during the synthe¬ 
sis. Abundances are also sensitive to the microturbulent 
velosity, so if the log g—vmicro expression used is not so 
accurate in this metallicity range, or the surface gravity 
is badly estimated, this propagates into systematically 
off abundances through v micro- 

The minimum [C/Fe] (Figure [6]) that can be measured 
from the CO lines strongly depends on temperature. Be¬ 
cause RGB stars in GCs generally have low carbon abun¬ 
dances, we can only set upper limits for a number of 
our stars. Our [C/Fe] values are an average of the de¬ 
rived abundances from five CO windows. As a result 
determining the upper limit is more challenging because 
CO lines in certain windows disappear faster with rising 
temperature than in others. We carefully checked ev¬ 
ery CO window fit and selected upper limits if the CO 
band head was not visible in more than three windows by 
looking at the flatness of the \ 2 fit around the minimum 
value. Because the derived abundance of N from the CN 
lines is anticorrelated with the value of [C/Fe] used, all 
stars with upper limits in [C/Fe] have also lower limits 
in [N/Fe]. These upper limits for [C/Fe] and lower limits 
of [N/Fe] are identified with open triangles in all figures. 

We see clear correlations between [C/Fe] and T e s 
in M13, M3, M2, and M5 when omitting upper lim¬ 
its. We interpret these as a sign of “deep mixing”, 
a nonconvective mixing process that causes steady de¬ 


pletion of surface carbon abundance and an enhance- 
rnent in ni troge n abundance in all low-mass RGB stars 
dGratton et alJ |2004[) . We derive only upper limits on 
[C/Fe] for M15, M92, M53 and NGC 5466 because the 
CO bands are quite weak at such low metallicity. With 
only upper limits, it is not possible to identify abun¬ 
dance trends in these clusters. However, previous studies 
(lAngelou et al.ll2012t iShetrone et ahll20 Kit [Mart ell et alJ 

2008) have found signs of the same deep mixing process 
in action in those clusters. 

The only cluster in which we see a sign of nitrogen 
enrichment with declining temperature is M13. Possible 
trends in nitrogen are not necessarily a result of astro¬ 
physics, because in our methodology small systematic 
errors in temperature can result in systematic errors in 
N abundance because deriving [N/Fe] is challenging from 
the CN lines. Other than the relatively large errors of 
N, one also has to be careful with possible correlations 
of C and N with temperature, because they can be gen¬ 
erated on the RGB as part of deep mixing and also in a 
previous stellar evolution event that are responsible for 
the pollution. 

In Figure [ 7 ] there is a weak trend, on par with the er¬ 
ror, visible in the Si abundance in M13, but not in the 
other clusters. While Ca does not show any correlation 
with temperature, it does show temperature-dependent 
scatter in most of the GCs. The three Ca lines used 
in our analysis are generally weak, and get significantly 
weaker above 4700 K, which leads to higher errors re- 
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Figure 7. [Si/Fe], [Ca/Fe], and [Ti/Fe] as a function of photometric T e fj for all ten clusters. The error bars represent our final combined 
uncertainties from Table 4. Trends in Ti for M2 and M5 are removed using the plotted lines when calculating the internal scatter (see 
Section 5.1 for discussion). 

lated to increased sensitivity to the uncertainty in the 
continuum placement. Titanium, unexpectedly, shows 
a decline with decreasing temperature in M2, M5 and 
(marginally) in M107 in Figure [ 7 ] We suspect that in¬ 
accuracies in our analysis of the hotter star spectra are 
driving this apparent trend, since the S/N is lower for 
those stars than for the cooler, brighter giants. Because 
these trends in Mg (Figure 0, Si and Ti (Figure [7J only 
show up in a handful of clusters, we believe that they are 
a result of difficulties in data analysis for certain lines in 
certain stars and not any systematic mishandling in, e.g., 

Vmicro estimation. 

We choose to fit these various observed trends with 
a linear equation, regardless of their origin, in order to 
explore the scatter around the trend. We fit lines to 
[Mg/Fe] in M107 and M71, [C/Fe] in M13, M2, M3 and 
M5, and [Ti/Fe] in M2 and M5. We are only using these 
fits to determine the internal scatter, which we define as 
the standard deviation about the fitted line. Through¬ 
out the rest of this paper we use the abundance values 
directly, except when discussing the CN anticorrelations 
in Section 6.4. 

5.2. Scatter and Errors 

An internal abundance scatter significantly larger than 
the estimated errors suggests an astrophysical origin. In 
fact, significant scatter in N, O , Mg a nd A1 is well doc¬ 
umented in literature dGratton et al.ll20l2l) . The large 
star-to-star variations in these elements are the result of 
the CNO, Ne-Na, and Mg-Al cycles, which play an im¬ 






portant role in the nuclear fusion processes at the early 
stages of RGB (CNO) and AGB (Ne-Na, Mg-Al) evo¬ 
lution. Figure [5] shows the internal scatter calculated 
as the standard deviation around the mean values (in 
some cases around a fitted linear equation, as explained 
earlier), and the final combined estimated uncertainties 
from Table 4, as a function of cluster average metallic- 
ity. We note that the small number of stars analysed in 
NGC 5466 limits our ability to make a detailed analysis 
of C, N, and O in this cluster: because we were only able 
to measure the abundance of these elements in the three 
stars that have temperatures below 4500K, the scatter is 
unrealistically higher than the calculated uncertainty. 

The star-to-star scatter in [Fe/H] is quite similar to 
our measurement uncertainties, indicating that no signif¬ 
icant metallicity variations in the clusters are detected in 
this study. Recently, lYong et all (|2014T ) have discovered 
three distinctive groups with different iron abundances 
in M2. The first group has a dominant peak at [Fe/H]=- 
1.7, while the second and third have smaller peaks at 
1.5 and 1.0, but the membership for the latter group is 
not conclusive. We see no such behaviour in M2. The 
discrepancy between these t wo stu d ies is most likely the 
result of a selection effect. lYong et all (i2014jl se lected 
stars that belong to a second RGB ( Lardo et al.ll2012j ) 
within the cluster, while ou r sam pl e select ion were based 
upon previous observations (iZasowski et alJ[2(Tl3h and all 
of our stars belong to the main RGB of this cluster. 

The a—elements Si and Ti also show constant dis- 


[Si/Fe], [Ca/Fe], [Ti/Fe] as a function of T eff , T eff < 5400K 
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Figure 8. 

dots). 


The internal scatter in each cluster from Table 6 (red dots) compared to our estimated final combined errors from Table 4 (blue 


tribution, but the scatter in [Ca/Fe] is larger than our 
measurement uncertainties for most clusters in our sam¬ 
ple. We attribute this to an increasing inaccuracy in 
[Ca/Fe] with rising temperature, as discussed in Section 
5.1, and not a real range in calcium abundance. Little to 
no range in iron and a—elements is what we expect for 
normal globular clusters, although the clusters in which 
these abunda nces d o vary are an intriguing puzzle (see, 
e.g., iMarino et al.l 120131 ). M107 stands out as having 
the largest discrepancy between the scatter of [Ca/Fe] 
and its estimated uncertainty. However, the large inter¬ 
nal scatter is not related to astrophysics, but to spectral 
contamination from telluric OH lines. Due to the clus¬ 
ter’s radial velocity, two of the three Ca lines lie on top of 
atmospheric OH lines. Because the telluric correction is 
not perfect, the Ca lines are compromised in a way that 
is not included in our uncertainty estimation. 

We do see scatter above the level of the uncertainties 
in some of the light-element abundances. In the case of 
carbon, we find that the scatter around the overall car¬ 
bon depletion trend for M13, M2, M3 and M5 (as shown 
by dashed lines in Figure [6]) is fairly small, although our 
measurement uncertainties get quite large at low metal- 
licity. We do not see clear car bon depletion i n Ml 07 or 
M71; however, as discussed in iMartell et al.1 (120081b the 
rate of carbon depletion due to deep mixing is lower in 
higher-metallicity stars. 

The range in [N/Fe] abundance (Figure© second panel 
from the top, in the left column) in M71, M107, M5, 
M3, M13, and M2 is clearly larger than the estimated 


uncertainties. However, the difference between the scat¬ 
ter and the uncertainty decreases with decreasing cluster 
metallicity, and below [Fe/H]=—1.7 (NGC 5466, M53, 
M92, and M15) our derived [N/Fe] values are all lower 
limits. The decreasing discrepancy is a result of larger 
uncertainties in more metal poor clusters, and the fact 
that the minimum N abundance of a cluster is slightly 
increasing as the average metallicity decreases (Figure 
©, while the maximum N is constant above [M/H]=- 
1.7. The uncertainties in M15 and M92 are significantly 
larger than in other clusters, and our measurements are 
not precise enough to judge the amount of N enrichment 
in these clusters. Fitting the CN lines in our defined 
CN window is also challenging. This is because they 
are usually blended with other molecular lines, and be¬ 
cause N abundance strongly depends on the derived C 
abundance, which is an upper limit in many cases. A 
more detailed analysis of the C-N anticorrelation and N 
bimodality can be found in Section 6.4. 

Our uncertainties of the [O/Fe] abundance (Figure© 
are fairly constant as a function of cluster metallicity. 
This is because the error is dominated by a strong sensi¬ 
tivity to the temperature used in the model atmosphere, 
and the internal scatter varies between clusters. Unfor¬ 
tunately, the Na lines available in the APOGEE spectra 
are too weak to investigate the O-Na anticorrelation in 
our sample, so we can only limit our discussion on the 
scatter of [O/Fe]. We use the later as an indication of 
the strength of O-Na anticorrelation in a cluster. In M71 
and M107 the O scatter is comparable to the uncertain- 
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Table 7 

Averages of Populations 


Cluster 

Pop. a 

N b 

R° 

[Fe/H] 

[C/Fe] 

[N/Fe] 

[O/Fe] 

[Mg/Fe] 

[Al/Fe] 

[Si/Fe] 

[Ca/Fe] 

[Ti/Fe] 

M15 

1 

10 

0.43 

-2.31 

-0.49 

0.78 

0.64 

0.20 

-0.19 

0.31 

0.08 

0.11 


2 

13 

0.57 

-2.26 

-0.33 

1.11 

0.43 

0.04 

0.75 

0.53 

0.24 

0.17 

M92 

1 

14 

0.30 

-2.23 

-0.42 

0.89 

0.67 

0.33 

-0.23 

0.41 

0.08 

0.09 


2 

33 

0.70 

-2.23 

-0.38 

1.03 

0.40 

0.05 

0.70 

0.47 

0.11 

0.09 

M53 

1 

8 

0.50 

-1.96 

-0.43 

0.91 

0.59 

0.38 

-0.11 

0.42 

0.23 

0.24 


2 

8 

0.50 

-1.95 

-0.54 

1.13 

0.55 

0.25 

0.84 

0.41 

0.23 

0.31 

N5466 

1 

6 


-1.83 

-0.56 

0.84 

0.64 

0.15 

-0.42 

0.28 

0.14 

0.29 


2 

2 


-1.79 

-0.56 

0.86 

0.62 

0.11 

0.31 

0.30 

-0.16 


M13 

1 

32 

0.40 

-1.50 

-0.44 

0.83 

0.56 

0.21 

0.04 

0.38 

0.24 

0.18 


2 

49 

0.60 

-1.49 

-0.58 

0.92 

0.12 

0.09 

0.99 

0.40 

0.28 

0.21 

M2 

1 

7 

0.39 

-1.46 

-0.45 

0.75 

0.56 

0.30 

-0.10 

0.35 

0.19 

0.26 


2 

11 

0.61 

-1.50 

-0.49 

1.01 

0.30 

0.24 

0.79 

0.35 

0.27 

0.27 

M3 

1 

39 

0.66 

-1.40 

-0.42 

0.59 

0.53 

0.18 

-0.08 

0.30 

0.13 

0.11 


2 

20 

0.34 

-1.38 

-0.54 

0.86 

0.16 

0.10 

0.79 

0.29 

0.10 

0.12 

M5 

1 

60 

0.49 

-1.24 

-0.41 

0.62 

0.41 

0.24 

0.06 

0.34 

0.20 

0.24 


2 

62 

0.51 

-1.24 

-0.50 

0.87 

0.16 

0.21 

0.64 

0.35 

0.20 

0.28 

M107 d 

1 

6 


-1.04 

-0.16 

0.47 

0.38 

0.36 

0.51 

0.50 

0.21 

0.12 


2 

6 


-1.00 

-0.25 

0.92 

0.29 

0.37 

0.53 

0.45 

0.22 

0.16 

M71 d 

1 

3 


-0.72 

0.00 

0.57 

0.53 

0.46 

0.44 

0.41 

0.24 

0.38 


2 

4 


-0.69 

-0.17 

1.16 

0.49 

0.44 

0.53 

0.38 

0.21 

0.45 


a Population 1 is denoted by red, while population 2 is denoted by blue in Figures 8—12. 
b N: the number of stars in each population. 

c R: the ratio of number of stars in a population and the overall number of stars analysed. 
d Populations are found using N instead of Al. 


ties, suggesting that the O-Na anticorrelation is not as 
extended as in the rest of the clusters. A clear spread 
is visible in M5, M3, M2, M13, M15 and M92, while 
in M53 the 0 abundances are constant and the scatter 
is significantly smaller than what is expected from the 
uncertainties. 

The effects of the Mg-Al cycle are very noticeable for 
all clusters except M71 and M107 (Figure [8l second pan¬ 
els from the bottom). The scatter in Mg abundance in 
those clusters is close to the estimated uncertainties af¬ 
ter taking the linear correlations with T e ff into account 
(Figure 0. In other clusters the scatter of [Mg/Fe] also 
closely follows the uncertainties, except for M15 and M92 
where the most Mg poor stars can be found. We see no 
Al enrichment in M107 and M71; however, as the aver¬ 
age metallicity decreases, the amount of internal scatter 
rapidly starts to deviate from the error, and we see a large 
spread in Al in all other GCs. A more detailed analysis 
of the Mg-Al anticorrelations and Al populations can be 
found in Section 6.2 and 6.3. 

Based on the examination of the spread of N and 
Al abundances, we conclude that all clusters possess 
multiple stellar populations either based on N or Al, 
or both. We examine these populations in more de¬ 
tail in the next section. The derived a—element abun¬ 
dances are also fairly constant in all globular clusters 
prese nted here, simil ar to what is reported in the litera¬ 
ture {Grafton et al.ll2004 ). 

6 . MULTIPLE POPULATIONS 

Most light elements show star-to-star variations in all 
GCs. These large variations are generally interpreted as 
the resul t of chemical feedback from an earlier ge neration 
of stars {Grafton et, al.ll200Tl : iCohen et al.l 12002ft . rather 
than inhomogeneities in the original stellar cloud these 
stars formed from. Thus, the current scenario of GC evo¬ 
lution generally assumes that more than one population 
of stars were formed in each cluster. The first generation 


of GC stars formed from gas that had been enriched by 
supernovae in the very early Universe, while the second 
formed by combining gas from the original star-forming 
cloud with ejecta from the first-generation stars. Only 
the fraction of first-generation stars contribute to the pol¬ 
lution, and the time scale of the formation of these second 
generation stars depend on the nature of polluters; it is a 
couple of hundred Myr in the case of intermediate-mass 
AGB stars, but it is only a few Myr for fast rotating 
massive stars and massive binaries. 

These first-generation polluters are thought to 
have been fast rotating massive stars dDecressin et al.l 
12007ft . or intermediate mass (M star > 3 Mq) 

AGB stars (iD’Ercole et all 12008ft . or massive binaries 
{de Mink et al.l 12009ft . Our data reveal the expected 
signatures of pollution from material enriched from the 
hot hydrogen burning cycles such as the CNO, NeNa, 
or MgAl cycles in all globular clusters in our sample. 

In this section we explore the various correlations be¬ 
tween these elements, and we discuss individual star for¬ 
mation, and/or pollution events suggested by separate 
groups found in the C-N and Mg-Al anticorrelations. 

6.1. Identifying Multiple Populations 

To separate the various populations und er stu dy, we 
follow an approach similar to that of jGratton et_al.1 
{2011ft . who used A'-means clustering {Steinhauslll956ft 
to identify multiple populations in to Cen. Here, we use 
t he ex treme-deconvolution (XD) method of iBovv et alJ 
{2011(1^1 to identify population groups and assign mem¬ 
bership. This method fits the distribution of a vector 
quantity, here the elemental abundances, as a sum of 
K Gaussian populations, whose amplitudes, centers and 
covariance matrices are left entirely free. The algorithm 
can be applied to noisy or incomplete data, thus making 

23 Code available at http: //github. com/jobovy/extreme-deconvolution . 





























16 


Table 8 

Statistics of Correlations 


Cluster 

A a 
[Mg/Fe] 

eA b 
[Mg/Fe] 

a c 

A 

[Si/Fe] 

eA 

[Si/Fe] 

CT 

M15 

0.16 

0.037 

4.3 

0.22 

0.055 

4.0 

M92 

0.28 

0.029 

9.7 

0.06 

0.038 

1.6 

M53 

0.13 

0.035 

3.7 

0.01 

0.055 

<1 

M13 

0.12 

0.014 

8.6 

0.02 

0.018 

1.1 

M2 

0.06 

0.038 

1.6 

0.00 

0.043 

<1 

M3 

0.08 

0.011 

7.3 

0.01 

0.024 

<1 

M5 

0.03 

0.008 

3.8 

0.01 

0.014 

<1 


a Difference of tile average Mg and Si abundance of populations. 
b The estimated error of the difference. 
c Detection significance in a. 


the best possible use of all available data. In the current 
application, XD’s main advantage is the latter, since we 
are not able to measure abundances for all nine elements 
in all stars in our sample. Similar to the A'-means, the 
number A of populations to fit is an input to the algo¬ 
rithm; XD itself does not determine the optimal number 
of components to fit the distribution. 

Briefly, the XD algorithm works by optimizing the like¬ 
lihood of the Gaussian mixture model of the data. The 
optimization proceeds by an iterative procedure consist¬ 
ing of repeated expectation (E) and maximization (M) 
steps. In the E step, each datum is probabilistically as¬ 
signed (a) membership in each population and (b) error- 
free values of each noisy or missing elemental abundance; 
in the M step, each Gaussian’s parameters are updated 
using its members’ mean abundances and covariance cal¬ 
culated from the error-free values obtained in the E step. 
These steps are repeated until the likelihood stops in¬ 
creasing to within a small tolerance. The algorithm is 
proven to increase the likelihood in each EM-step. 

We analyse each globular cluster in our sample us¬ 
ing this algorithm, with the [Mg/Fe],[Al/Fe], [Si/Fe], 
[Ca/Fe], [Ti/Fe] abundances as well as only [Mg/Fe] and 
[Al/Fe]. We find that the two separate analyses provide 
nearly identical results, indicating that Si, Ca, and Ti do 
not drive the populations. For M7I and MI07 we used 
[N/Fe] instead of [Al/Fe], because A1 does not show any 
spread in these two clusters. We run XD using uncertain¬ 
ties given by the ‘Final combined errors' section in Table 
4, and also without any uncertainties, and found that 
the two different runs produce identical results. We in¬ 
clude missing abundance measurements by using a large 
uncertainty for these measurements. Group membership 
for each star is determined using the best-fit Gaussian 
mixture by calculating the posterior probability for each 
star to be a member of each population based on its ele¬ 
mental abundances and their uncertainties; stars are then 
assigned to the population for which this probability is 
the largest. These membership assignments are typically 
unambiguous (probabilities > 99 % in most cases), with 
only a few cases for which the maximum probability is 
below 90 %. The abundance averages of populations are 
listed in Table 7. In the next few subsections we discuss 
the results from the population fitting in more detail for 
each individual cluster. 

6.2. The Mg-Al anticorrelation 


In order for the Mg-Al cycle to operate, high 
temp eratures above 70 million Kelvin are required 
dCharbonnel fe Prantzosl I2006H . Because current clus¬ 
ter main-sequence stars are unable to reach these tem¬ 
peratures, the high [Al/Fe] abundances we see in some 
globular cluster stars imply that a previous generation 
of higher-mass or evolved stars must have contributed to 
their chemical composition. 

Figure [9] shows [Mg/Fe] versus [Al/Fe] for all ten clus¬ 
ters in our study. The two XD populations are plotted in 
red and blue, and a representative error bar is included in 
the top right corner of each panel. As noted previously, 
stars in MI07 and M71 do not show a strong Mg-Al an¬ 
ticorrelation, and their XD population assignments are 
based on N abundances. The light blue points in the 
panels for M107 and M71 are stars that are warmer than 
4500 K, and thus have no CNO measurements. 

The extended distribution of Mg and A1 abundances 
apparent in Figure[!j]is typical for globular clusters, and it 
shows the influence of the Mg-Al fusion cycle, which con¬ 
verts Mg into Al. However, there is a variety in the struc¬ 
ture of this relationship that has not been thoroughly ex¬ 
plored before: in some clusters (M92, M53, NGC 5466, 
M2 and M3), there are two distinct abundance groups 
with a gap, while in other clusters (M15, M13 and M5) 
there is no gap. This result strongly suggests that there 
is diversity in the process of stellar chemical feedback and 
star formation in globular clusters, which may relate to 
the larger environment in which they formed. 

To date, the largest homogen eous study of M g and A l 
abundances was carried out by iCarretta et al.1 (l2009ail . 
We have several clusters in common, M15, M5, M107 and 
M7I which enables a direct comparison. M107 and M71 
do n ot show any a n ticorre lation in either of the stud¬ 
ies. iGarretta et all (12009a ) only have upper limits for 
Al and Mg in M15, while we were able to make direct 
measurements, but even with upper limits the two stud¬ 
ies show similar anticorrelations with largest spreads in 
both Al and Mg from the whole sample of clusters. The 
largest difference between the two studies can be seen in 
case of M5, where our values span a range_ o f —0.3 dex 
to +1.1 dex, while the Al abundances of IGarretta et al.l 
(l2009al) are limited to between —0.2 and 0.7 dex. The 
diffe renc e may be explained by a simple selection effect, 
as ICarretta et al.l (l2009al ) observed 13 stars in M5, while 
our sample size is 122 and covers almost the full extent 
of the RGB. 

We discuss the extent of anticorrelation through dif¬ 
ferences in the average of the Mg abundance between 
the populations, divided by the estimated uncertainty in 
this average Mg difference. We use the final combined 
errors from Table 4. Table 8 lists the errors of the av¬ 
erage differences with values of how much a detection is 
each difference. We find that in M15 (4.3<r), M92 (9.7cr), 
M13 (8.6cr), and M3 (7.3cr) the differences are statis¬ 
tically significant and the Mg-Al anticorrelation exists. 
Two clusters, M53 (3.7cr) and M5 (3.8tr), also show large 
a detections, however, we would like to use a more con¬ 
servative approach to the detection of anticorrelations 
out of a concern that our errors may be underestimated, 
given the degeneracies between stellar parameters and 
abundances in our analysis. M2 stands out as having no 
statistically significant Mg-Al anticorrelation. 

The summed abundance A(Mg+Al) is expected to be 
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Mg-AI anti-correlations, T eff < 5400K 
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Figure 9. Mg-AI anticorrelations. Colors mark the populations found with the XD code, first population is denoted by red, and the 
second-generation denoted by blue points. Because [N/Fe] was used to identify populations in M107 a nd M 71, not all stars could be 
catalogued, these are denoted by light blue. The same colors denote the same stellar groupings in Figures |91131 


A(Mg+AI) as a function of T pff 
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Figure 10. The combined abundance of A(Mg+Al) as a function 
of effective temperature. 

constant as a function of T e g- when material is completely 
processed through the Mg-AI cycle, and that is what our 
results show in Figure fll)l The only exception is M107, 
but the slight correlation is due Mg correlating with tem¬ 
perature (see Figure 0, and this trend was not removed 
here. 

A1 is expected to correlate with elements enhanced by 
proton-capture reactions (N, Na; Figure [TT| and anti¬ 
correlate with those depleted in H-burning at high tem¬ 
perature (O, Mg; Figure fllTl . The Al-0 anticorrelations 
and Mg-0 correlations can be clearly seen in our data 


in Figure [TTl for all clusters except M53, M107 and M71. 
The small number of stars observed in NGC 5466 does 
not allow us to investigate the correlations in that clus¬ 
ter in detail. However, the slight correlation of A1 with 
Si that can be seen in Ml5 in Fig. |TT| is the evidence of 
28 Si leaking from the Mg-AI cycle, which is an intriguing 
result. Si participating in the light-element ab undance 
patter n was first rep orted in NGC 6752 by lYong et al.1 
(|2005H . Since then, iCarretta et al.l (|2009af) have also 
found Si enhancement correlated with A1 in NGC 2808. 
Those authors’ interpretation was that it is only in low- 
metallilcity clusters, where the AGB stars burn slightly 
hotter, or in high-mass clusters, where the chemical en¬ 
richment is more efficient, that a Si enhancement will be 
observed in second-generation GC stars. The difference 
in Si abundance between the two XD-identified popu¬ 
lations, listed in in Table 8, is significant relative to the 
error on that measurement. The Si-Al correlation in M15 
is also accompanied by a Si-Mg anticorrelation (Figure 
EED, which is further evidence of 28 Si bein g prod u ced by 
hot bottom burning in AGB stars ([Karakas fc Lattanziol 

nm). 

6.3. The spread of Al abundances 

The spread of Al abundances (Figure 0 also increases 
significantly below [Fe/H]=—1.1, thus we can conclude 
that high mass, low metallicity AGB polluters that are 
able to process material at high, 60—70 million K, play an 
increasingly larger role in more metal-poor cluster. This 
behaviour of Al abundances was previously observed by 


























18 


LL 


OJ 

LL 


<D 


OJ 


OJ 

LL 

O) 


OJ 

LL 

O) 


0> 

LL 

< 


OJ 

LL 


OJ 


<D 


OJ 

LL 

O) 


<D 

LL 

D) 


1 

0 

1 

0 

1 

0 

1 
0 
1 

0 

1 

0 

0 10 10 10 10 1 


1 

0 

1 

0 

1 

0 

1 

0 

1 

0 

1 

0 


[Si/Fe] [Si/Fe] [Si/Fe] [Si/Fe] [Si/Fe] 



Si and O (anti-)correlations 

t -1-n=pi ^~r -1-m———r 


“ i i n 

. ^ . 

- M 15 V* 

- N 7078 1 ,. 

%• 

1 1 T 

- ■ 

- M 92 

- N 6341 ^ 

w ' NGC 5466 '. 

• 

• 

.N. 

- * f- 

- M 53 * 

- N 5024 *,* 

- M 13 

- N 6205 A* 

- m 

- M 2 \ 

- M 3 T 


m 

s€ 


- N 7089 . 

•> 

- N 5272 0 . 

V 1 

- M 5 *N 5 ^ 

- M 107 , N 6171 

- M 71 . N 6838 

-. J . 

•• 

1 1 1 

. • • 

V 

* 

I_ 1 _ 1 _ 

i 1 ' 

y 

I 1 1 

-II 1 1 - 

— 

in 

y 

- M 15 , N 7078 

_l_ 1 _ 1 _L_ 

- M 92 . N 6341 

_l_ 1 ___ 1 _L_ 

- NGC 5466 

- M 53 , N 5024 

j_ 1 _ i_ l 

- Ml 3 . N 6205 

i \ i i 

— 1 - 1 - 1 

s 

f • 

_ 1 _ 1 _I 

1 1 I 

• 

_ 1 _ 1 _I 

I 1 1 

s 

I_ 1 _ 1 _ 

I 1 1 

-*C W . 

i 

I i i 

- :• - 

- M 2 , N 7089 

- M 3 , N 5272 

- M 5 , N 59*04 

- M 107 , N 6171 

- M 71 , N 6838 

l-SH 


m 


m 

. <r . 

*• • *• . 

- Ml 5 , N 707 & * - 

: 

- M 92 , N 63 “ 

A* 

- NGC 5466 

- M 53 , N 5024 

- M 13 , N^f 



\m 

\=\ 



m 


** 

•Sr 

- M 2 , N 7089 

_1 _ 1 _ 1 _L 

- M 3 , N 5272 

J_ 1 _ 1 _L 

- M 5 , N 5904 

J_ 1 _ 1 _L 

- M 107 , N 6171 

J_ 1 _ 1 _L 

-M 71 .N 6838 


0 10 10 10 10 1 
[O/Fe] [O/Fe] [O/Fe] [O/Fe] [O/Fe] 


Figure 11. Si and O (anti)correlations with Al, N, and Mg. For explanation of color coding, please see description of Figure [9] 

iCarretta et all (l2009af l. but the larger range of Al values diet this behavior ( Ventura et alif200ll 12013 1. The high- 
presented in this paper make this correlation clearer. mass AGB stars reach higher temperatures at the bottom 

Theoretical AGB nucleosynthesis modeling indeed pre- of the convective envelope; i.e., stronger HBB and ad- 
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CN anti-correlations, T eff < 4520K 
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Figure 12. CN anticorrelations. Correlations of [C/Fe] with temperature associated with deep mixing were removed in clusters marked 
by ***. Upper limits are denoted by open triangles. 


vanced (Mg-Al) nucleosynthesis occurs with decreasing 
metallicity. We are seeing very advanced Mg-Al nucle¬ 
osynthesis in the most metal poor clusters such as M15 
and M92, while the most metal-rich clusters like M107 
and M71 do not show this, as theoretically expected, if 
the high-mass AGB stars are the polluters. Also, this is 
corroborated by the Al-0 anticorrelation; the most Al- 
rich stars are O-poor showing the effects of very strong 
hot bottom burning (HBB), because HBB proceeds com¬ 
pletely and destroys O. 

An other related issue is that HBB is activated for 
lower masses with decreasing metallicity. High-mass and 
very low-metallicity AGB stars don not exists in GCs to¬ 
day due to their extremely short lifetimes, but this trend 
seems to be confirmed both theoretically (jVentura et al.1 
12001112013T ) and observationally at least for solar metal- 
licities down to that of the Magellanic Clouds. From 
the observational point of view, we know that high-mass 
AGB stars in the Small Magellanic Cloud ([Fe/H]=—0.7) 
activate HBB for progenitor masses (M >3 M 0 ) lower 
th an their solar met allici ty co unterparts (M > 4 M 0 ) 
(jGarcia-Hernandez et al.l 120061 l2009h . Thus, at the low¬ 
est metallicities in M15 and M92 we also would expect 
more HBB AGB stars (i.e., with several degrees of HBB 
and Mg-Al nucleosynthesis), because the minimum stel¬ 
lar mass to activate the HBB process (and advanced Mg- 
Al nucleosynthesis) decreases with decreasing metallicity. 
Thus, we conclude that more polluted material would be 
present at the lowest metallicities. 

Because of these reasons, we believe that our results 


support and add some evidence to the high-mass AGBs 
as GCs polluters. 

6.4. C-N 

Carbon, nitrogen and sometimes oxygen are influenced 
by two independent processes in GC giants: a primordial 
anticorrelation, with the same first-generation sources as 
the O-Na and Mg-Al patterns; and also stellar evolution, 
driven by circulation between the hydrogen-burni ng she ll 
and t he surface (|Sweigart fe Mengellll979t lAngelou et al.l 
120131 . As discussed in the previous section, we see large 
star-to-star variations in N abundance in all clusters; 
however, our inability to measure C abundances below 
[Fe/H]=-1.7 dex leads us to only having upper limits, 
and therefore we restrict the discussion of C-N anticor¬ 
relations to the more metal-rich clusters. 

Figure fl2l shows [N/Fe] versus [C/Fe] for all clusters. 
Deep mixing has to be taken into account when dis¬ 
cussing the CN anticorrelations, so in M13, M2, M3, and 
M5, where we see deep mixing clearly, the correlation of 
[C/Fe] with T e g was removed. Because deep mixing is 
not visible in M107 and M71, we used the original values. 
A clear anticorrelation cannot be seen in any of the clus¬ 
ters. While M107 and M71 show weak anticorrelations, 
neither of these is statistically significant. The anticorre¬ 
lation itself may be obscured by the relatively large errors 
for [C/Fe] and [N/Fe]. As mentioned in Section 5.1, we 
see correlations of [C/Fe] with temperature in M13, M3, 
M2, and M5, but these are not accompanied by increas¬ 
ing [N/Fe] at the same time. We believe that this is due 
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Figure 13. Upper panels: The sum of C, N and O as a function of 
effective temperature. Upper limits are denoted by open triangles. 
For explanation of color coding, please see description of Figure 
8 . Lower panels: The sum of C, N and O as a function of [Al/H], 
A clear correlation is visible in M13, M2, M3, and M5, for more 
discussion see Section 6.4. 

to our inability to measure C and N accurately: most of 
our abundance determinations are upper limits or close 
to them. If the anticorrelations exist, they must span 
ranges smaller than our uncertainties, which are always 
higher than 0.14 dex in [C/Fe] and 0.12 dex in [N/Fe]. 
It is interesting to note that in M3, M107 and M71 there 
is a clear CN weak and CN strong group of stars, while 
in all more metal poor clusters, N abundances fill out a 
broad distribution. 


6.5. C+N+O 

We have C, N, and O abundances available for a 
large number of stars, thus we are able to investi¬ 
gate the C+N+O content in each cluster in detail. 
Studies from the literature showed that the C+N+O 

content in glo bular cluste rs _ are fairly constant to 

within 0.3 dex divans et all 119991 iCarretta et al.l 120051 


Yongj+jd 


Yong et al 


S et al.l 1199a) . except for N1851 where 
find a spread of 0.57 dex. According to 
the large spread in C+N+O in N1851 is probably 
attributed to larger than usual pollution from lower mass 
AGB stars tha n in other c lust ers, b ut these results were 
questioned by IVillanova et all (I2010H as they did not see 
spread of C+N+O larger than 0.3 dex. 


We also see near constant C+N+O in our sample, 
which is consistent with the material in these stars having 
undergone CNO cycling in the first-generation of stars in 
the RGB phase. In Figure [Rj] we show A(C+N+0) as a 
function of T e ff. The spread is significantly larger than 
what is reported in literature, between 0.4 and 0.6 dex 
for M13, M2, M3, and M5. The spread in these clusters 
is at the level of what has been found in N1851, but we 
think this is mostly associated with large uncertainties of 
[N/Fe] measurements. We do not find clear separation in 
the amount of A(C+N+0) between the two populations 
found in any of the clusters. 

As previously mentioned, there are currently two lead¬ 
ing models to explain the nature of the polluters: the 
first assumes that intermediate-mass AGB stars that are 
thermally pulsing and undergoing hot bottom burning 
expel mate rial to the intra-clu ster medium by strong 
mass loss (IVentura et all 120011 1 , while the second as¬ 
sumes th at the po llution come s from hot, fast rotat¬ 
ing stars dDecressin et al.l 1200711 . According to the first 
theory, intermediate-mass HBB-AGB stars could pro¬ 
duce o r no t pro duce large C+N+O variations (see, e.g., 
lYong et al.ll2009t and references therein). This is because 
the C+N+O predictions (which basically depend on 
the number of third dredge-up episodes) for these stars 
are extremely dependent on the theoretical modelling, 
the convective model a nd the ma ss lo ss prescriptio n 
used (see, e.g,. IDAntona fe Ventural120081 : IKarakas et all 
1201211 . A (C+N+0)-Al correlation is not clearly visible 
in the lower panel of Figure 10. 

According to the first theory, AGB stars may explain 
a large spread in A(C+N+0) because they are expected 
to produce a substantial increase in the C+N+O abun¬ 
dance as t hey enhance Na and Al and deplete O and Mg 
( Yong et all 12009 1. This should result in a (C+N+O)- 
A1 correlation, which is not clearly visible in the lower 
panel of Figure fl3l Weak correlations in M107 and M71, 
and weak anticorrelations in M13 and M3 are visible, but 
they are probably the result of large uncertainties in C 
and N. More accurate measurements of [N/Fe] will help 
to clarify the possible (C+N+0)-Al correlation. 

7. SUMMARY 

We investigated the abundances of nine elements for 
428 stars in ten globular clusters using APOGEE DR10 
spectra. A homogeneous analysis of these GCs has not 
been accomplished previously, something that APOGEE 
is uniquely able to do because it can observe all the bright 
GCs in the northern hemisphere. A semi automated code 
called autosynth was developed to provide abundances 
independent of those derived by ASPCAP. Our main goal 
was to examine the stellar populations in each cluster 
using a homogeneous dataset. Based on our abundances, 
we find the following: 

1. From the examination of the star-to-star scatter in 
a—element abundance, only O and Mg show an abun¬ 
dance scatter larger than our measurement errors. How¬ 
ever, this is expected, because O and Mg are part of 
the well-known globular cluster light-element abundance 
anticorrelations. 

2. Population analysis using [N/Fe] for M107 and 
M71 and [Al/Fe] in the rest of the clusters confirms that 
each cluster can be divided into two discrete populations, 
though in M107 and M71 that division is clearest in the 
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N-Al plane, whereas in the other clusters it is best made 
in the Mg-Al plane. 

3. The anticorrelation of Mg and A1 is clearly shown 
in M15, M92, M13 and M3, and to a lesser extent in 
M53 and M5. Interestingly, M2, which is also a metal 
poor cluster, does not have a Mg-Al anticorrelation. The 
increased number of stars in our sample compared to 
the literature enables us to discover more Al-rich stars, 
making the anticorrelation clearer in this dataset than in 
previous studies. 

4. The spread of A1 abundances increases with decreas¬ 
ing cluster metallicity. This is in agreement with t heoret¬ 
ical AGB n ucleosynthesis predictions by lYentura et al.1 
(|2001l I2013D . This suggests that high-mass HBB AGB 
stars, which are able to enrich A1 while destroying O, are 
more common polluters in metal-poor clusters than in 
metal-rich ones. 

5. A correlation between A1 and Si in M15 indi¬ 
cates particularly high-temperature hydrogen burning 
in the stars that contributed the abundance pattern of 
the second-generation stars. This is consistent with Si 
enhancement in globular cluster stars as discussed in 
the literature, where it appears to be confined to low- 
metallicity and high-mass clusters. 

6. Besides accessing northern GCs, APOGEE is also 
unique among other spectroscopic surveys, because it 
is capable of measuring abundances of CNO. The to¬ 
tal abundance of C+N+O is found to be constant in 
each cluster, within our errors. Our data set do not 
show an unambiguous correlation between A(C+N+0) 
and [Al/H]. In principle, this finding is not against the 
idea of intermediate-mass HBB-AGB stars as the pri¬ 
mary source of chemical self-enrichment in globular clus¬ 
ters. However, more precise measurements of C and N 
abundances will clarify this issue, and we expect that 
ASPCAP will be able to provide these in the future. 
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